Knudsen layer formation in laser induced thermal desorption 
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Laser induced thermal desorption of Xe atoms into vacuum from a metal surface following the nano-second 
pulsed laser heating was investigated by the time-of- flight (TOF) measurement. The desorption flow was 
studied at a wide range of desorption flux by varying the initially prepared Xe coverage O (1 ML = 4.5 x 
10 18 atoms/m 2 ). At 8 = 0.3 ML, the TOF of Xe was well represented by a Maxwell-Boltzmann velocity 
distribution, which is in good agreement with thermal desorption followed by collision-free flow. At O > 0.3 
ML, the peak positions of the TOF spectra were shifted towards the smaller values and became constant at 
large O, which were well fitted with a shifted Maxwell-Boltzmann velocity distribution with a temperature Tb 
and a stream velocity u. With Td fixed at 165 K, u was found to increase from 80 to 125 m/ s with increasing 
O from 1.2 to 4 ML. At O > 4 ML, the value of u become constant at 125 m/ s. The converging feature of 
u was found to be consistent with analytical predictions and simulated results based on the Knudsen layer 
formation theory. We found that the Knudsen layer formation in laser desorption is completed at Knudsen 
number Kn < 0.39. 

PACS numbers: 47.45.-n, 47.70.Nd, 68.43.Nr, 68.43.Vx 



I. INTRODUCTION 

Laser induced desorption of atoms and molecules 
from solid surfaces is a vital phenomenon to investi- 
gate fundamentals such as surface electronic structures 
and dynamics^— Besides, laser desorption is an essential 
technique for pulsed laser deposition used in thin film 
growth-r— and mass spectrometry of protein employ- 
ing matrix-assisted laser desorption ionization. 15 When 
desorption flux is small, the velocity distribution of des- 
orbed atoms is directly governed by the desorption mech- 
anism. When the desorption flux is large enough, on the 
other hand, the post-desorption collision between des- 
orbed particles may become significant and modify the 
velocity distribution in the vicinity of the surface after 
the desorption. 

Manifestations of the collision effect in laser induced 
desorption have been reported both by experiments and 
simulations as the modifications of the angular and ve- 
locity distribution of desorbing atoms and molecules 
Cowin et al. investigated the angular dependence of the 
translational temperature of D2 desorbed from tungsten 
surfaces under a pulsed laser irradiation.— The trans- 
lational temperature of D2 desorbed in the surface nor- 
mal direction was higher than those in oblique directions. 
They attributed the variation of the translational tem- 
perature to the collision effect. Noorbatcha et al. used 
the direct Monte Carlo simulation of desorbing atoms 
to investigate the collision effect. They showed that 
even in the sub-monolayer regime the collision notice- 
ably modifies the final angular, velocity and rotational- 
energy distributions.— ~— However, there has not been 



a ) E-mail: a-ikeda@iis.u-tokyo.ac.jp 

b ) E-mail: fukutaniOiis . u-tokyo. ac. j p 



any model that can quantitatively estimate the degree 
of modification by the post-desorption collision in laser 
desorption so far. 

Knudsen layer formation theory has been developed 
in rarefied gas dynamics to model the steady flow of 
the strong evaporation from the surface*^— As shown 
in Fig. [I] the initial velocity distribution of thermally 
desorbed species at the surface is well described by a 
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FIG. 1. (a) A schematic sketch of the Knudsen layer formed 
above the surface following the laser induced thermal des- 
orption. The desorption flux is assumed to be large enough 
for the post-desorption collisions to take place and to further 
form the Knudsen layer. The arrows denote the velocity vec- 
tors of desorbed atoms at the surface and at the end of the 
Knudsen layer, (b) Schematic plots of the velocity distribu- 
tion of the desorbed atoms in the z (surface normal) direction 
at the surface and at the end of the Knudsen layer. Ts , Tk 
and u denote the surface temperature, the translational tem- 
perature of the desorbed atoms at the end of the Knudsen 
layer and the stream velocity, respectively. 
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"half-range" Maxwell-Boltzmann velocity distribution. 18 
In the Knudsen layer theory, as a result of intensive post- 
desorption collisions, the half-range velocity distribution 
at the surface is thermally equilibrated to a full-range 
Maxwell-Boltzmann velocity distribution with a stream 
velocity at some distance from the surface. This thermal- 
ization layer is defined as the Knudsen layer. The theory 
analytically predicts for monoatomic gas that the ratio 
of the translational temperature at the end of Knudsen 
layer Tk to the surface temperature Tg and the Mach 
number of the desorption flow at the end of the Knudsen 
layer become 0.65 and 1.0, respectively 21 ' 22 



Substrate on a cold head 




Kelly and Dreyfus discussed that the Knudsen layer 
formation theory may be applicable to the pulsed des- 
orption flow with a large desorption flux. However, it 
is not straightforward because the pulsed desorption in- 
volves complex time evolution of the density and velocity 
distributions of desorbed atoms; 23 Sibold and Urbassek 
have shown by means of the Monte Carlo simulation of 
the Boltzmann equation that the pulsed desorption flow 
at an intense flux is well characterized by the above val- 
ues predicted by the Knudsen layer formation theory. 
Although previous experimental studies recognized the 
collision effect in laser desorptionj^ir— the formation of 
the Knudsen layer in laser desorption has been discussed 
only by theory and simulations i 21 i 22 ' 29 So far, any exper- 
imental confirmation of Knudsen layer formation in laser 
desorption has not been presented. For an experimen- 
tal verification of the theory, a systematic observation of 
the translational temperature and stream velocity of the 
desorption flow as a function of the desorption flux is 
strongly required. 

In the present paper, we investigated the laser induced 
thermal desorption (LITD) of Xe from an Au(001) sur- 
face by means of the time-of-flight (TOF) measurement 
as a function of the wide range of desorption flux by 
varying the initial Xe coverage 9. We found that at O 
close to ML, the TOF was well analyzed by a Maxwell- 
Boltzmann velocity distribution. Hence, the desorption 
at O close to ML is rationalized by the thermal des- 
orption followed by the collision-free flow. At close to 
monolayer, we observed that the peak positions of the 
TOF spectra shift towards smaller values. Assuring that 
the desorption is only thermally activated, we regard this 
modification of the TOF as the manifestation of the col- 
lision effect. At larger 0, the peak positions of the TOF 
spectra become constant. We deduced the Mach number 
M of the desorption flow to be 0.96 at large under the 
assumption that Tk/Tg = 0.65. The obtained value of 
M and the saturating behavior of u at > 4 ML both 
well coincide with the previously reported values by the 
Knudsen layer theory and simulations. The facts sug- 
gest the formation of the Knudsen layer in LITD at large 
0. Furthermore, we tentatively estimated the Knudsen 
number of the initial desorption flow. 



FIG. 2. A schematic drawing of the experimental setup for the 
time-of-flight measurement of Xe following the laser induced 
thermal desorption from Au surfaces. C. Amp., PD and BS 
denote fast current amplifier, photo diode and beam splitter, 
respectively. 



II. EXPERIMENT 

The Au specimen and the physisorbed Xe layers were 
prepared in the following manner. An Au disk with 001 
orientation was cut out from a single-crystal rod and was 
chemically and mechanically polished. The Au disk was 
put into an ultra high vacuum (UHV) chamber with a 
base pressure of 2 x 10~ 8 Pa. The sample surface was 
cleaned by several cycles of Ar + sputtering and anneal- 
ing at 700 K in the UHV by the electron bombardment. 
The cleanliness of the specimen was confirmed by the 
observation of a clear reconstructed (5 x 20) pattern 3 ^ by 
low-energy electron diffraction and no contamination in 
the Auger electron spectrum. The sample temperature 
was monitored with a Chromel-Alumel type thermocou- 
ple directly spot-welded to the side of the Au disk. The 
disk was cooled down and kept at 23 K by a closed-cycle 
He compressor type refrigerator during all LITD experi- 
ments. 

The Xe gas was dosed onto the Au surface by back- 
filling the chamber. The Xe coverage was monitored 
by the LITD yield. The LITD yield showed a linear de- 
pendence on the Xe exposure in the present experimental 
condition. The saturation coverage at the sample tem- 
perature of 80 K was defined to be 1 ML = 4.5 x 10 18 
atoms/m 2 as in Ref. l3ll was varied by varying the 
dosage of Xe gas onto the Au surface. We note that 
the desorption flux is defined as devided by the mean 
desorption period r. On the basis of the result employ- 
ing the LITD simulation, t was shown to be constant at 
about 4 ns at concerned in the present study; 32 Hence, 
the desorption flux can be varied by varying 0. 

The TOF measurement of the LITD of Xe from an Au 
surface was carried out in the following manner. The ex- 
perimental setup is schematically shown in Fig. [5J An 
ArF excimer laser (Lambda Physik) was used as a pulsed 
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FIG. 4. The stream velocity u of the desorbed Xe atoms 
from Au surfaces following the pulsed laser irradiations as a 
function of the Xe coverage 0. The data are obtained by an- 
alyzing the time-of-flight spectrum with the shifted Maxwell- 
Boltzmann velocity distribution. The solid line is a guide for 
eyes, u become constant at above 4 ML which is denoted by 
the arrow. 



recorded using the oscilloscope with a single acquisition 
mode. 



FIG. 3. Time-of-flight (TOF) spectra of Xe from Au surfaces 
following pulsed laser irradiations. The Xe coverage are 
10, 6.7, 1.7 and 0.3 ML from top to bottom, respectively. 
The dashed curves and the solid curves are the Maxwell- 
Boltzmann (M-B) velocity distribution and the shifted M-B 
velocity distribution fitted to the experimental results, respec- 
tively. The vertical dashed line and the arrows indicate the 
peak position of each spectrum. 



laser source. The wavelength was 193 nm and the time 
duration was 8 ns. The incidence angle was 25° from 
the surface normal direction. The irradiated area on the 
sample disk was ~0.1 cm . The laser pulse was guided 
on to the surface through an aperture and a quartz win- 
dow. The laser power absorbed by the sample was fixed 
at 80 mJ/cm 2 in all experiments taking the reflectivity 
into consideration. At this laser power absorbed by the 
sample, the LITD of Xe is the dominant desorption pro- 
cess as shown in Ref. [32|. The reflected light escapes 
the UHV chamber through another quartz window. A 
quadrupole mass spectrometer (QMS: Balzers QMA125) 
was placed in the surface normal direction at a distance 
of 10 cm from the Au specimen. The QMS was tuned 
for a high-sensitive detection of Xe with a low mass se- 
lectivity. The ion current of the QMS was amplified with 
a fast current amplifier (Keithley 427) and sent to an os- 
cilloscope (Tektronixs: TDS620B). The desorption of Xe 
from the Au(001) surface was induced with only one laser 
pulse in each LITD experiment. The respective TOF was 



III. RESULTS 

We obtained a series of TOF spectra of desorbed Xe 
from Au surfaces following pulsed laser irradiation as 
shown in Fig. [3] The TOF spectrum of Xe from the 
Au surface at O = 0.3 ML is shown at the bottom of Fig. 
[3] The vertical dashed line at 0.42 ms shows the peak 
position of the spectrum. The spectrum was well ana- 
lyzed with a Maxwell-Boltzmann velocity distribution in 
a flux from a thermal source described as^ 

(TThl) \ 

where m, k and v are the mass of a Xe atom, the Boltz- 
mann constant and the velocity of Xe atoms, respec- 
tively, and A and Td are fitting parameters. In the 
analysis, we convert Eq. ([l) to a TOF function for 
the density sensitive detector, which results in the form 
f(t)dt = at- 4 exp(~bt- 2 )dt. By fitting Eq. {TJ to the 
TOF of Xe at 6 = 0.3 ML, we obtained T D of 255 K as 
illustrated by the dashed curve at the bottom of Fig. [3l 
At the top and middle of Fig. |3l the TOF spectra 
of desorbed Xe at O = 10 to 1.7 ML are shown. We 
note that the peak positions of those TOF spectra, which 
are indicated by the arrows, are shifted towards smaller 
values with increasing O compared with that of 6 = 0.3 
ML. It was found that the peak positions of the TOF 
spectra become constant at about 0.34 ms at O > 5.0 
ML. 
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IV. DISCUSSION 

A. Desorption mechanism 

We first discuss that the observed desorption proceeds 
only via thermal activation. In the previous study^ we 
investigated the same system by varying the photon flu- 
ence at 9 = 1 ML. At the small fluence region where 
the laser power is below 20 mJ/cm 2 , we observed a des- 
orption of Xe only at a photon energy of 6.4 eV and no 
desorption at 2.3 eV. We assigned this photodesorption 
as a non-thermal desorption via a transient Xe - forma- 
tion. At a large fluence region where laser power exceeds 
60 mJ/cm 2 we observed, on the other hand, that thermal 
desorption is the dominant desorption process regardless 
of the photon energy. These observations are confirmed 
by comparing the experimentally observed Tb and Xe 
desorption yield with the results of LITD simulations. 
Therefore, we regard that the initial desorption is only 
thermally activated in the present experiment. 



B. Analysis of TOF spectra 

In order to analyze the TOF spectra at 9 > 0.3 ML, 
we introduce a shifted Maxwell-Boltzmann velocity dis- 
tribution which is characterized by u. Generally, the ve- 
locity distribution in a desorption flux shows an angular 
distribution and is expressed by the elliptical distribu- 
tion which involves an angular dependent translational 
temperature as described in Ref. |3J. In the present ex- 
periment, the detector was fixed in the surface normal 
direction. Therefore, the translational temperature only 
in this direction is detected in the present experiment. 
In this case, the shifted Maxwell-Boltzmann velocity dis- 
tribution in a flux from a thermal source is expressed 

W-^expj 2kJh\ () 

Under the condition of u = 0, Eq. ([2]) is identical to Eq. 
([T]). In the analysis, Eq. (J2J was also converted to a TOF 
function for the density sensitive detector. 

The TOF spectra at 9 > 0.3 ML were well analyzed 
using Eq. <j2j) - However, we found it difficult to obtain 
a unique set of Td and u. The spectra were well repro- 
duced with an arbitrary value of Td from 165 to 310 K by 
adjusting u from 125 to m/s. We note that, with a fixed 
value of Td between 165 and 310 K, the obtained value of 
u monotonically increases with increasing O from 0.3 to 
4 ML and becomes constant at 9 > 4 ML. For example, 
Fig. |4] shows u as a function of 9 obtained by using Eq. 
@ with Td fixed at 165 K, of which the fitting curves 
are shown as solid curves in Fig. [3J 



C. Moderate desorption at 6 ~ ML 

At 9 ~ 0, the observed feature of the TOF can be 
rationalized with a model employing thermal desorption 
followed by collision-free flow. At the instance of thermal 
desorption, the desorbed Xe gas is in thermal equilibrium 
with the surface at the temperature Tg. Here, it is rea- 
sonable to consider the half-range Maxwell-Boltzmann 
velocity distribution for desorbed Xe expressed as2i 

(mv 2 \ 

v z > (3) 

where A and v are a normalization factor and the velocity 
vector of Xe atoms in the Cartesian coordinate, respec- 
tively. We note that in Eq. ((3]) the velocity component 
in the z direction has a distribution only at v z > Opi 
assuming the z axis to be the surface normal direction. 
Given the desorption period is much shorter than the 
flight time and the irradiation diameter is much smaller 
than the flight distance, Eq. ^ can be, as a velocity dis- 
tribution in a flux from a thermal source, transformed to 
the same form as Eq. ([TJ. Thus, we notice that Td = Tg 
at 9 ~ 

For a qualitative analysis, we estimated the surface 
temperature during the present LITD experiment using 
the first order desorption rate equation, assuming that 
the desorption activation energy is 240 meV— and that 
the desorption period is 4 ns i 35 ' 36 The calculated results 
show that the required surface temperature is 260 K, 
which is in good agreement with the obtained value of 
T s = 255 K at 9 = 0.3 ML. The agreement indicates 
that the model of the thermal desorption followed by 
the collision-free flow well describes the observed TOF at 
9 ~ 0, in good accordance with the results by Cowin 14 
and Wedler^. 



D. Intensive desorption at 6 > 4 ML 

At 9 > 4 ML, we discuss that the desorption flow is de- 
scribed as an intense flow, where the post-desorption col- 
lision modifies the velocity distributionj 15 ' 16 Several the- 
oretical studies have speculated that the so-called Knud- 
sen layer is formed in the vicinity of the surface as a result 
of collisions! 21 ' 22 ! 29 ' 37 As schematically drawn in Fig. [T] 
in the Knudsen layer model, the initial velocity distri- 
bution is in thermal equilibrium with the surface at Tg 
except that there is no distribution at v z < as in Eq. 
([3]) . As a result of a significant number of post-desorption 
collisions, the half-range velocity distribution at z = be- 
comes thermally equilibrated at some distance from the 
surface to a full-range Maxwell-Boltzmann velocity dis- 
tribution with a stream velocity u i 18 ' 19 ' 21 The model also 
requires that, at the end of the Knudsen layer, the tem- 
perature of the flow becomes identical (Tk) in all degrees 
of freedom. Hence, at the end of the Knudsen layer, the 
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TABLE I. Mach number M obtained in the present study 
and previously reported values by simulation and theory on 
the Knudsen layer. 



Method 


M 


Condition 


Experiment^ 
Simulation 29 
Simulation 38 
Theory^^ 


0.96 
1.0 
1.0 
0.99 


Pulsed flow 
Pulsed flow 
Steady flow 
Steady flow 



a Present study 

velocity distribution may be described a a 18 ' 19 i 21 
/(v)e?v = Aexp 



m { v l + v l + (Vz ~ u k) 2 } 



2kT K 



dv, (4) 



where Vi, uk and Ik are the velocity component in the 
Cartesian coordinate, stream velocity and the tempera- 
ture at the end of the Knudsen layer, respectively. Un- 
der the condition that the flight distance is sufficiently 
longer than the thickness of the Knudsen layer besides 
the conditions described above, Eq. Q is converted, as 
a velocity distribution from a thermal source, to the iden- 
tical form to Eq. @. Thus, we see that Td = Ik and 
u = uk at > 4 ML. Equation ([2]) was shown to well 
describe the experimental results of the TOFs at > 4 
ML. This fact, along with the convergence feature of u, 
indicates that the experimental result at > 4 ML may 
be rationalized with the Knudsen layer formation model. 



E. Mach number of the flow 

We further examine the results at > 4 ML by 
quantitatively comparing the obtained results of «k with 
previous theoretical reports on Knudsen layer forma- 
tion in the steady strong evaporatio n 18 ' 19 and the pulsed 
desorption- 29 - from plane surfaces. We compare uk in 
terms of the Mach number M at the end of the Knudsen 
layer. M is defined as 



UK 



jkT K ? 



(5) 



where c and 7 are the local velocity of sound and the 
heat capacity ratio of the gas, respectively, the latter of 
which is 5/3 for monoatomic gas as in the present case. 

Ytrehusi 8 - and Cercignanii^ formulated M and Tk/Tg 
by finding a solution to the Boltzmann equation under 
the Knudsen layer formation model assuming the con- 
servation of the particle number, momentum and the en- 
ergy flux between the surface and the end of the Knudsen 
layer. Sibold and Urbassek estimated M and Xk/Ts us- 
ing the Monte Carlo simulation of the Boltzmann equa- 
tion in one-dimension for both pulsed flows 9 , and steady 
fiow^ 8 . conditions. Those reports have stated that at the 
end of the Knudsen layer, M should always be close to 
unity and that Xk/Ts is at around 0.65. 



As discussed in IV. B, it was difficult to unambiguously 
determine the value of Tjy. Here, we assume the relation 
of Tk/Ts = 0.65 following the theoretical studies and 
evaluate the Mach number. With the value of Ts = 255 
K, we obtain the value of T K to be 165 K. With T K = 165 
K and Eq. (|5]l. we obtain the value of M to be 0.96. As 
can be seen in Table I, the previously reported values 
of M show a quantitative agreement with the obtained 
value of M in the present study. This indicates that the 
observed trend of u at > 4 ML is consistent with the 
Knudsen layer formation theory. 21 

With adiabatic expansion, M should well exceed 
unity,2L22d£ We note that we observed the saturation 
of M at around unity at > 4 ML in the present ex- 
perimental conditions. This indicates that the adiabatic 
expansion is practically absent, although the slight in- 
crease of u at > 4 ML may be due to the adiabatic 
expansion. 



F. Knudsen number of the flow 

In order to generalize the present result, we tentatively 
consider mean gas density n, mean free path A and Knud- 
sen number Kn in the vicinity of the surface at the mo- 
ment of desorption as a function of 0. Kn is defined 
by^ 



Kn=^ = —f}- 

H y/2nal z 



(6) 



where l z and a are the representative length of the system 
and the Van der Waals collision cross section of Xe (1.0 x 
10~ 19 m 2 ), respectively. 

Here, we simply regard that l z is the mean thickness 
of the gas cloud above the surface. It, then, reads that 
h — v z T , where v z is the mean thermal velocity of des- 
orbed Xe in the z direction and r is the mean desorption 
period. We further take v z = y / 2kTs/Trm and estimate 
that n = Q/l z . By substituting n in Eq. (|6]), we obtain 
a simple relation Kn = a/(y/2Q) as introduced in Ref. 
I29I . and a are the relative coverage = 0/©s with 
the monolayer saturation coverage ©s and the ratio of 
the area occupied by an atom at ©s to a described as 
a = 1/(<j©s), respectively. 

Now, we see that Kn depends simply on © and a. 
Therefore, we can discuss the required condition for the 
formation of Knudsen layer in LITD for general systems 
in terms of Kn. We fixed t = 4 ns. On the basis of the 
result employing the LITD simulation, r showed little 
dependence on in the region concerned in the present 
study. 

In Table HH we summarized the obtained n, A and Kn 
as a function of 0. Kn can also be understood as an 
inverse of mean collision times per each atoms. In the 
present study, we observed the formation of the Knudsen 
layer at > 4 ML, which corresponds to Kn < 0.39 and 
to more than ~ 2.6 collisions per desorbing atom. The 
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TABLE II. Mean gas density n, mean free path A and Knud- 
sen number Kn in the vicinity of the surface at the moment 
of laser desorption as a function of G. 



e (ml) 


n (10 24 atoms/m 3 ) 


A (nm) 


Kn 


0.3 


3.3 


2114 


5.2 


4.0 


45 


159 


0.39 


10 


111 


63 


0.16 



manifestation of collision effects is observed at > 0.3 
ML, which corresponds to Kn > 5.2. The results coin- 
cide with the previous theoretical work by Sibold and 
Urbassek 2 - that the Knudsen layer is formed at = 2.5 
ML and that the collision effect already appears while 
the Knudsen layer is not formed at = 0.25 ML, respec- 
tively. The obtained values of Kn also agree well with the 
previous experimental observation of the collision effect 
in LITD of D 2 /W at = 1 ML reported by Cowin et 
alJ^. They also coincide with the theoretically estimated 
collision number of 2.9 per atom at = 1 ML reported 
by Noorbatcha et al.^S- 

Lastly, we note the observed feature at ~ 1 ML, 
where the Knudsen layer is not formed. Thus, the veloc- 
ity distributions are not necessarily well fitted with Eq. 
@. Nevertheless, they were very well analyzed with Eq. 
([2]) in all at ~ 1 ML. The result is in good agreement 
with the simulated results by Sibold and Urbassek^ al- 
though any analytical formulation for the flow at ~ 1 
ML has not been presented so far. Hence, further theoret- 
ical and experimental works are required for elucidating 
the mechanism for the modification of velocity distribu- 
tion at ~ 1 ML due to the post-desorption collision as 
pointed out in Ref. [29j. 



V. CONCLUSION 



In conclusion, we observed the TOFs of Xe from a 
Au(001) surface by LITD at a wide range of 0. At 
~ 0.3 ML, the TOF is well analyzed by the Maxwell- 
Boltzmann velocity distribution at 255 K indicating ther- 
mal desorption followed by the collision-free flow. At 
> 0.3 ML, the TOF is well analyzed by the shifted 
Maxwell-Boltzmann velocity distribution, which is re- 
garded as the manifestations of the post-desorption col- 
lision effects. At > 4 ML, the value of u under the as- 
sumption that Tk/Ts = 0.65 became constant at around 
125 m/s, which corresponds to M — 0.96, being in good 
agreement with the prediction by the Knudsen layer for- 
mation theory. In the LITD experiment, we found that 
the collision effect already appears at Kn ~ 5.2 and that 
the Knudsen layer is formed at Kn < 0.39, which cor- 
responds to a mean collision number of greater than 2.6 
per atom for general systems. 
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